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RESUMEN 

Estudiamos las propiedades de una muestra de 27 Nucleos Activos de Galaxias 
(NAGs) con perfil de doble pico en las lineas anchas de emision de baja ionizacion, 
tornados del SDSS (DR4) con un valor medio de onu B ~ 3002 ± 139km -s -1 . 
Nuestro primer resultado es que todo el espectro de emision en el rango de 4100A 
a 5800A se puede ajustar muy bien mediante un modelo de emision en un disco de 
acrecion eliptico, suponiendo que las lineas de H/3, Fell, H7 y HeIIA4686A tienen 
el mismo perfil de doble pico que las lineas de Ha. Un analisis similar se ha hecho 
en un articulo previo para el objeto SDSS J2125-0813. Los mejores ajustes indican 
que la emision optica de las lineas de Fell en los emisores dbp (emisores de lineas 
de doble pico, por sus siglas en ingles) se originan en la misma region del disco que 
las lineas de Balmer. Encontramos que algunas correlaciones conocidas para los 
NAGs "normales" se cumplen tambien para los dbp. Sin embargo, este resultado 
debe tomarse con cautela porque el numero de objetos es pequeo y porque tenemos 
un sesgo de seleccion hacia objetos con emision de Fell intensa. Mostramos que 
para los dbp, la masa del agujero negro central parece tener mayor influencia en las 
propiedades de la emision de Fell que la taza de acrecion. Tambien encontramos 
que todos los dbp de nustra muestra son cuasares radio callados, con excepcion 
de ocho objetos con R r > 1 de acuerdo a la definition de Ivezic et al. (2002) y 6 
objetos no descubiertos por FIRST. 

ABSTRACT 

We study the Fell properties of double-peaked broad low-ionization emission 
line AGN (dbp emitters) using a sample of 27 dbp emitters from SDSS (DR4), with 
mean value <JHa B ~ 3002±139km • s _1 . Our first result is that the line spectra in the 
wavelength range from 4100A to 5800A can be best fitted by an elliptical accretion 
disk model, assuming the same double-peaked line profiles for H/3, Fell, H7 and 
HcIIA4686A as that of double-peaked broad Ha for all the 27 dbp emitters, except 
the object SDSS J2125-0813 which we have discussed in a previous paper. The best 
fitted results indicate that the optical Fell emission lines of dbp emitters originate 
from the same region in the accretion disk where the double-peaked Balmer emission 
lines originate. Some correlations between Fell emission lines and the other broad 
emission lines for normal AGN can be confirmed for dbp emitters. However, these 
results should be taken with caution due to the small number of objects and the 
bias in selecting strong Fell emitters. We show that for dbp emitters, BH masses 
seems to have more influence on Fell properties than dimensionless accretion rate. 
We also find that the dbp emitters in the sample are all radio quiet quasars except 
one dbp emitter with R r > 1 according to the definition by Ivezic et al. (2002) and 
6 objects undiscovered by FIRST. 
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1. INTRODUCTION 

The optical Fell emission lines are one of the 
most important properties of some Active Galac- 
tic Nuclei (AGN). Because of the large amounts of 
Fell emission lines (half-filled 3d-shell leads to thou- 
sands of transitions for Fell emission lines) (Moore 
& Merill 1968; Netzer 1988; Wills et al. 1980; Wills, 
Netzer & Wills 1980), it is difficult to clarify many 
properties of these lines, such as their loci of emis- 
sion, the mechanism of their excitation, the corre- 
lation between Fell lines and other emission lines. 
The bright object I Zw I (PG 0050+124) has very 
strong Fell emission lines which can be accurately 
measured. And thus it has been used to build Fell 
templates for the subtraction of these lines when re- 
ducing quasar spectra in order to properly fit the 
continuum (Philips 1978; Oke & Laucr 1979; Boro- 
son & Green 1992; Laor et al. 1997; Marziani et al. 
2003a). 

A detailed study of line-related correlations be- 
tween Fell and other properties of AGN can be found 
in Sulcntic, Marziani & Dultzin-Hacyan (2000): 
there is an outstanding strong anticorrelation be- 
tween EW(FeII4570)/EW(H/3) and FWHM(if/3) in 
the context of the so called Eigenvector 1 (El) pa- 
rameter space. These authors have defined a popu- 
lation A, FWHM(H/3 broad )< 4000km • s" 1 , and pop- 
ulation B, FWHM(H/3f, road )> 4000km • s"\ which is 
a cleaner distinction for AGN with strong and weak 
Fell emission lines than Radio loud objects vs. Radio 
quiet objects. A preliminary definition of Eigenvec- 
tor 1 can be found in (Boroson & Green 1992). An 
updated and more complete version of it includes the 
interpretation of the observed trends in this parame- 
ter space in the context of their physical drivers: BH 
masses and accretion rate (Sulentic et al. 2006). 

There are two main kinds of models which can 
reproduce the observed shape and equivalent width 
of Fell emission lines: Photoionization models with 
microturbulence (Korisra et al. 1997; Bottorff et 
al. 2000; Netzer 1985; Baldwin et al. 2004) in the 
context of the so called LOC model (Baldwin 1995) 
and collisionally excited models (Grandi 1981, 1982; 
Kwan et al. 1995; Dumont et al. 1998; Baldwin 
et al. 2004). A simple photoionization scheme is 
not adequate to reproduce the observed Fell lines 
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(Bergeron & Kunth 1980; Collin-Souffrin et al. 1986; 
Joly 1987; Collin-Souffrin, Hameury & Joly 1988) 
because of the enhanced Fell emission due to the 
needed high density (> 10 n cm~ 3 ) of the emitting 
clouds (Ferland & Person 1989). Thus, the accre- 
tion disk and/or the region near the center produced 
by shocks along the radio jet provide a natural high 
density environment for Fell emission lines. Sigut 
& Pradhan (2003) have examined this assumption 
for the typical physical conditions of AGN. In their 
paper, the following excitation mechanisms for Fell 
emission lines have been considered: continuum fluo- 
rescence (Phillips 1978, 1979), collisional excitation 
(Joly 1991), self- fluorescence among the Fell tran- 
sitions and fluorescent excitation by Lya and Ly/3 
(Pcnston 1988; Sigut & Pradhan 1998; Verner et al. 
1999). 

The most commonly accepted outcome from the 
study of Fell emission lines is that the line width of 
these lines is the same as that of the other broad 
low-ionization emission lines such as H/3, except 
for several special objects such as HE 1249-0648 
and HE 1258-0823 (FWHM( J ff/?)>>FWHM(FeII)) 
(Marziani et al. 2003a, 2003b). The best way to de- 
termine where do the Fell emission lines originate, 
the accretion disk or the base of the jet, is to study 
the properties of Fell lines from a special kind of 
AGN which emit double-peaked broad low-ionization 
Balmer emission lines (dbp emitters). These double- 
peaked broad lines are believed to originate from 
the accretion disk near the central Black Hole (BH) 
(Storchi-Bergmann, Nemmen, et al. 2003; Storchi- 
Bergmann, Eracleous et al. 1997; Storchi-Bergmann, 
Eracleous & Halpern 1995; Eracleous et al. 1997; 
Storchi-Bergmann, Baldwin et al. 1993; Chen et al. 
1989a, 1989b, 1997; Halpern et al. 1996; Antonucci 
et al. 1996; Sulentic et al. 1990; Shapovalova et 
al. 2001; Gilbert et al. 1999; Hartnoll & Blackman 
2002; Karas et al. 2001). Other models have been 
considered to interpret the origination of doubloe- 
peaked broad lines, such as the Binary Black Hole 
model (Begelman et al. 1980; Gaskell 1983) and 
bipolar outflowing model (Zheng et al. 1990), but 
they have been proved unsuccessful to explain most 
dbp emitters. 

In a previous paper we have reported a 
flat-spectrum radio quiet quasar, SDSS J2125- 
0813, which has a broad double-peaked H/3 ( 
FWHM(i?/3)~ 15000km- s- 1 ), and strong Fell 
emission lines at optical bands that have exactly 
the same line profile as that of the broad double- 
peaked H/3 line (Zhang, Dultzin-Hacyan & Wang 
2006a, hereafter paper I). In this paper, we select 
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a whole sample of AGN with both double-peaked 
broad Balmer lines and strong optical Fell emission 
lines to study the Fell properties of dbp emitters. 
There are two famous samples of dbp emitters, one 
consists of 23 objects which are nearly all LINERs 
from radio galaxies (Eracleous & Halpern 1994, 2003; 
Eracleous et al. 1995), and the other sample is made 
up of 112 objects of which 12% arc classified as LIN- 
ERs (Strateva et al. 2003) selected from SDSS DR2 
(York et al. 2000). Recently, we have built a larger 
sample of dbp emitters from SDSS DR4 (Zhang et 
al. 2006b), which includes more than three hundred 
dbp emitters. Because of the convenience to get the 
spectra of SDSS, we select dbp emitters with appar- 
ent Fell emission lines from the two samples in SDSS 
(Strateva et al. 2003; Zhang et al. 2006b). In section 
2, we list the data of our sample. Section 3 presents 
the results and then the discussions and conclusions 
follow in Section 4. The cosmological parameters 
H Q = 75 km s^Mpc" 1 , n x = 0.7 and fl m = 0.3 
have been adopted here. 

2. THE SAMPLE 

There are 112 dbp emitters selected from SDSS 
DR2 by Strateva et al. (2003) and more dbp emit- 
ters are being selected from SDSS DR4 by Zhang et 
al. (2006b). Because we select dbp emitters from 
the main catalogs of galaxies and of quasars classi- 
fied by pipeline of SDSS, with different criteria from 
those used by Strateva et al. (2003), and because 
in this paper we simultaneously select dbp emitters 
with apparent and strong Fell emission lines, even 
when we use DR4, our sample is smaller (some dbp 
emitters selected by Strateva et al. (2003) are not 
included in our sample). 

The detailed selction criteria will be described in 
a forthcoming paper (Zhang et al. 2006b). We give 
here a description of the steps we followed: First, 
we subtracted the stellar components from the ob- 
served spectra by the PCA template method (Hao 
et al. 2005; Li et al. 2005), if necessary. Then the 
emission lines, especially Ha+[NII]A6548, 6583A, are 
fitted by gaussian functions. At least two gaussian 
functions, one broad and one narrow, are applied 
for each permitted emission line, and one gaussian 
function for each forbidden emission line. According 
to the best fitted results by Levenberg-Marquardt 
least-squares minimization technique, the AGN with 
broad emission lines are selected according to the fol- 
lowing criteria: the line width (er), and flux of the 
broad component of Ha should be at least 5 times 
larger than the measured errors, and also the line 
width g of the broad component of Ha should be 



larger than 800km • s . The value of the summed 
squared residuals divided by degrees of freedom, x\ > 
can also be obtained. The value \\ near to 1 repre- 
sents whether the parameters of the model are sig- 
nificant for the fit to the emission lines. Because of 
the complex line profiles of double-peaked emission 
line, one broad gaussian function is not the better 
choice to fit the broad component, the values of \\ 
for objects with double-peaked broad Ha should be 
larger than 1. Thus multi-gaussian functions: Four 
broad gaussian functions, are applied for broad emis- 
sion lines again. Then We obtain the other value of 
x\ for each AGN with broad emission lines. For 
objects with double-peaked broad Ha, the value of 
x\ should be much closer to 1. We then select our 
dbp emitter candidates according to the following 
criterion: \i ~ X 2 > 1) which indicates the observed 
line profiles can not be properly fitted by one broad 
gaussian function. Finally, we check the dbp emitter 
candidates one by one by eye, according to whether 
there is apparent double peaks. We reject the candi- 
dates which have asymmetric line profiles of Ha or 
have multi-gaussian components with the same cen- 
ter wavelength (the separation of peaks is less than 
lOA). 

According to the observed spectra of SDSS, we 
can easily select the dbp emitters for which the Fell 
emission lines are evident, especially in the wave- 
length range from 4100l to 5800A. We find 27 dbp 
emitters with apparent Fell emission lines including 
the object SDSS J2125-0813, 12 dbp emitters are 
selected from the sample of Strateva et al. (2003) 
and 15 new dbp emitters are selected from SDSS 
DR4. We show the spectra in Fig 1 (the spectra of 
SDSS J2125-0813 can be found in Paper I, here we 
do not show it again). According to the properties 
of the spectra, all the 27 dbp emitters are quasars. 
The other 26 objects have both double-peaked Ha 
and H/3 emission lines (because we have a limit on 
the redshift of dbp emitters selected from SDSS less 
than 0.33). This allows us a more convenient way 
to measure and estimate the Fell emission lines as 
discussed below. The objects are listed in Table 1. 

3. RESULTS 

3.1. Measurement of the line parameters of 
emission lines 

Due to the presence of double-peaked broad Ha, 
in order to get the best fitted results for narrow 
emission lines, there are three to four broad gaus- 
sian functions for the broad components of Ha, 
one single narrow gaussian function for each nar- 
row emission line in the wavelength range between 



4 



ZHANG, DULTZIN-HACYAN & WANG 



6250A and 6850A and a power law for the contin- 
uum. Thus, there are three to four broad gaussian 
functions for broad Ha and seven narrow gaussian 
functions for [NII]A6548, 65831, [SII]A6719, 6732A, 
[OI]A6300, 6363A and narrow component of Ha. 
Here, we limit the center wavelength of each doublet 
to the same redshift, and limit the velocity disper- 
sion of gaussian function of each doublet to the same 
value in velocity space. Furthermore, we fix the flux 
ratio of [NII]A6548l to [NII]A6583A to 0.33, limit 
the line width of narrow emission line (au ne ) less 
than 400km • s _1 and limit the line width of broad 
emission lines (aline) larger than 500km • s _1 . The 
last best fitted results are shown in Fig 1. 

According to the best fitted results for narrow 
emission lines, we can get the broad components of 
Ha in the wavelength range from 6250A and 6850A. 
This provides a better way to measure the Fell emis- 
sion lines at optical bands and a better way to es- 
timate the model parameters of the accretion disk 
model. Because of strong Fell emission lines, it is 
much more difficult to obtain the complete line pro- 
file of H/3. However, we can estimate broad double- 
peaked H/3 by the line profile of Ha after the sub- 
traction of narrow emission lines. Before proceed- 
ing further, we should first subtract the continuum 
near H/3. We subtract the continuum according to 
the best fitted results for the points near 4100A and 
near 5800A. The continuum near H/3 is also shown 
in Fig 1. 

Once we obtain the spectra in the wavelength 
range between 4100A and 5800A after the subtrac- 
tion of the continuum, we can fit the spectra as we 
have done in paper I, under the assumption that Fell 
emission lines come from the same regions as broad 
Balmer emission lines: 

fx =k-Hfi x Ha + k Hl x Ha + k FeII FeII+ 

k Fell I ^ FeIH + k MgI ^ MgI + H Pnarrow + 
H"/ n arrow + [Oil I] 

(1) 

where ku ne is the flux ratio of the emission line 
to broad Ha. The narrow component of H/3, H7 
and [OIII] doublet can be fitted by narrow gaussian 
functions. If necessary, there is a broad component 
with the same line profile as that of broad Ha for 
HcIIA4686A, and two extended gaussian functions 
for the extend wings of the [OIII] doublet. The best 
fitted results are shown in Fig 2. The parameters for 
the flux ratio of different emission lines are listed in 
Table 2. 

From the fitted results shown in Fig 2, we notice 



that there are some objects for which the emission 
lines near H/3, especially near 4600A, can not be bet- 
ter fitted. The main reasons are perhaps the follow- 
ing: First, we do not consider the contributions of 
narrow emission lines of H7+[OIII] and CVI emis- 
sion lines, since they are weak emission lines. The 
second reason the accuracy of the continuum sub- 
traction. We fit the continuum under H/3 by one 
power law function. However, From the study of 
composite spectra of quasars (Francis et al. 1991; 
Zheng et al. 1997; Berk et al. 2001), especially the 
study of composite spectra of AGN selected from 
SDSS, we can see that the continuum should be best 
fitted by two power laws with a break of <~ 5000A 
Thus, applying only one power-law to fit the contin- 
uum under H/3 should lead to some uncertainty. 

After the subtraction of narrow emission lines 
near Ha, we can obtain the parameters of the el- 
liptical accretion disk model (Eracleous et al. 1995) 
by the Levenberg-Marquardt technique and adjust- 
ing the 8 free parameters (MPFIT package in IDL). 
There are eight free parameters in the elliptical ac- 
cretion disk model, the inner radius ri„, the outer 
radius r out , the inclination angle of the accretion 
disk i, the local broadening velocity dispersion cr, 
the slope of the line emissivity q, the eccentricity e, 
the orientation angle of the accretion disk </> and an 
amplitude factor k. The best fitted results are shown 
in Fig 3. The parameters of the accretion disk model 
for each dbp emitter are listed in Table 2. 

We select elliptical accretion disk model rather 
than circular disk plus spiral arm model, due to the 
following considerations. For elliptical disk model, 
there are mainly four parameters to determine: The 
observed line profile, ri n , r out , i and e, however in 
the circular disk plus spiral arm model, there are an- 
other more four parameters that dominate the ob- 
served line profile. To some extend, we think the 
uniqueness of model parameters for elliptical accre- 
tion disk model is better than circular disk plus spiral 
arm model. Also, we do not select a simple circu- 
lar accretion disk model, because the model can not 
be applied for dbp emitters which have brighter red 
than blue peak. 

Last, we measure the line parameters of broad 
Ha, such as the line width of Ha, UHa B ■ The line 
parameters of the broad component of Ha can be 
measured by: 



A 



JXfxdX 
IfxdX 

J X 2 fxd\ 
IfxdX 



J(X- Xp) 2 f\dX 
IfxdX 



(2a) 
(2b) 



THE PROPERTIES OF OPTICAL FEII EMISSION LINES OF DBP EMITTERS 



5 



We find that the mean value of f\ after the subtrac- 
tion of the continuum is zero for the nearby points on 
the continuum, however, the mean value (A — Ao) 2 /a 
for those points further out is not zero, and becomes 
larger for the points selected from more extend dis- 
tance from the center of the emission line, unless we 
use a smoother background, after the subtraction of 
the continuum, for which the value of each point is 
equal to zero. Thus, we should select cautiously the 
wavelength to measure the second moment according 
to the equation above. The wavelength range is se- 
lected as follows: after the subtraction of the contin- 
uum and the bad pixels, we select the last point for 
which the flux is equal to the minimum value (some- 
times, the minimum value of the flux is not equal to 
zero) of the points on the blue side of the emission 
line and select the first point for which the flux is 
equal to the minimum value of the points on the red 
side of the emission line. Due to the composite na- 
ture of the line profile of double-peaked emission, we 
think it is more accurate to measure the second mo- 
ment of the broad Ha using the equation above, than 
to measure the value of FWHM of Ha according to 
the standard definition of FWHM, because the flux 
density of the broad component is the largest one 
in the composite emission line after the subtraction 
of the narrow emission lines. The value of onu B is 
listed in Table 1. 

3.2. Properties of Fell emission lines of dbp 
emitters 

First, we check the correlation between line width 
of broad Balmer emission lines and strength of Fell 
emission lines EW(FeII)/EW(H/3). We show the cor- 
relation in Fig 4. The mean value of the second 
moment of broad Ha is about 3002 ± 139km ■ s _1 . 
As pointed out above, we use the second moment 
of broad Ha rather than FWHM(Ha). The 
Spearman Rank correlation coefficient is -0.39 with 
Pnuii ~ 5%. Here, the data include the object 
SDSS J2125-0813, for which the second moment of 
broad Ha is estimated from the broad component of 
H/3 reproduced by the accretion disk model. The 
results indicate that the anti-correlation between 
EW(FeII)/EW(H/3) and line width of broad lines is 
also valid for dbp emitters. 

The correlation between Fell properties and the 
line width of Balmer broad lines reflects, to some 
extent, effects of the line-of-sight. In the case of 
the correlation between the Fell properties and the 
width of the narrow lines, which are considered as 
tracers of stellar velocity dispersion, the effect of the 
bulge-BH mass on the properties of Fell emission 



may be more neatly deconvolved. A detailed descrip- 
tion about this is shown in the next subsection. Fig 
5 shows the correlation between EW(FeII)/EW(H/3) 
and line width of narrow lines. Here, we select the 
line width of the [OIII] emission line. As described 
above, there two gaussian components are needed for 
[OIII] emission line: one normal (core) component 
and an extended one. We use the core components 
as described in Greene & Ho (2005a). The Spear- 
man Rank of the correlation coefficient is -0.52 with 

Pnull ~ 0.6%. 

Furthermore, we check the radio properties of 
dbp emitters with apparent Fell emission lines. We 
find that there are 8 dbp emitters in our sample tar- 
geted and 13 dbp emitters covered by FIRST/NVSS. 
Thus, we can get the radio flux density or the upper 
limit value of radio flux density at 20cm. According 
to the definition of radio loudness, by Ivezic et 
al. (2002): 

Ri = log(F 2Qcm /F l ) = 0.4 x (m, - t) (3) 

where m, is one of the SDSS magnitudes and 
t is 20cm AB radio magnitude defined as f = 
—2.5 • log(F 2 ocm/3631Jy), we can calculate ra- 
dio loudness for each dbp emitter in our sample. 
There are only one dbp emitter with R r > 1, 
SDSS J075407.95+431610.5, which is one radio loud 
quasar. The other 20 dbp emitters are radio quiet 
quasars. Thus dbp emitters with stronger Fell emis- 
sion lines are apt to be found in radio quiet quasars, 
which in agreement with that of El analysis of single 
peaked BRL AGN (e. g. Sulentic et al. 2003). 

It is also convenient to check the correlation be- 
tween the continuum luminosity and the luminosity 
of Ha which is confirmed by Greene & Ho (2005b) 
using a sample of AGN selected from SDSS, be- 
cause in our code, the flux of Ha is one of funda- 
mental parameters to determine the flux of other 
emission lines, such as the flux of double-peaked 
H/3, x fifa, the flux of Fell emission lines, 

KfcII x fnu- One of the main considerations is that 
the correlation reflects some physical parameters of 
BLRs of AGN, such as the covering factor of BLRs, 
the proportional contribution of the ionization en- 
ergy to the emission of broad emission lines etc.. The 
correlation is shown in Fig 7. The Spearman Rank 
correlation coefficient is ~ 1 with P nu u ~ 0. The 
correlation for dbp emitters obeys the same relation 
as that for normal AGN, Ln a ~ i 1157 o. Here, 

5100A 

the luminosity of Ha also includes the contributions 
of [Nil] doublets, besides the contributions of broad 
and narrow components of Ha, because sometimes it 
is difficult to separate the narrow component of Ha 
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from the components of [Nil] for some dbp emitters. 
Because of the smaller contributions of [Nil] emis- 
sion lines, it is not affect the result between Ln a 
and L ° . 

510CL4 

In the case of dbp we can get, from the best fitted 
results of the accretion disk model, the inclination 
angle of the accretion disk. The plot of the strength 
of Fell vs. the inclination angle of the accretion disk 
is shown in Fig 8. The Spearman Rank correla- 
tion coefficient is 0.43 with P nu ii ~ 6%. We selected 
only 20 objects with accurate inclination angles i, 
i > 1.5 x i error . The inclination angle derived from 
the accretion disk model has a somewhat large error, 
specially when the angle is near 90 degrees. We can- 
not confirm a trend for all objects, because, if we only 
consider the objects with inclination angle less than 
60 degrees, the coefficient is 0.16 with P nu u ~ 63%. 
On the other hand, if the trend between the strength 
of Fell emission lines and the inclination angle shown 
in Fig. 7 was true, it would indicate that the objects 
with broader Balmer lines have stronger Fell emis- 
sion lines, which is in contradiction with the much 
stronger correlation between strength of Fell emis- 
sion lines and the line width of broad Ha. And in 
contradiction with all previous results found from El 
parameter space analysis, e.g. Sulentic, Marziani & 
Dultzin-Hacyan (2000). Thus, the trend most be a 
fake one. This is due not only to the fact that values 
for the inclination have large errors for large angles, 
but also due to the fact that inclination alone cannot 
govern the Fell strength, it is always convolved with 
BH mass and accretion rate, as discussed below (see 
also Marziani et. al. 2003a) 

3.3. Influence of Black Hole mass and Accretion 
Rate 

The most accurate way to estimate the central 
BH masses is by the stellar velocity dispersion as in 
Gebhardt et al. (2000); Ferrarese & Merritt (2001); 
Tremaine et al. (2002). For quasars, it is not pos- 
sible to determine the BH masses by stellar veloc- 
ity dispersion. Another convenient way is under the 
assumption of Virialization (Greene & Ho 2005b; 
Ovcharov et al. 2005; Wu et al. 2004; McLure & 
Jarvis 2004; 2002; Marziani et al. 2006; Onken et al. 
2004; Peterson et al. 2004). However, there are var- 
ious caveats for the measurement of the full width 
at half maximum (FHHM) of broad low-ionization 
emission lines, such H/3, especially for dbp emitters. 
A detailed discussion of the line components that do 
reverberate can be found in Sulentic et al. (2006). 
Actually, these authors prefer to use the Fell lines 
(assuming they are virialized as we prove here to 



determine BH masses). Unfortunately, for double- 
peaked broad emission lines, it is not reasonable to 
fit the line by one gaussian function or to measure the 
value of FWHM for double-peaked emission lines. 

However, there is a strong correlation between 
stellar velocity dispersion and line width of low- 
ionization narrow emission lines, which has been 
studied by many authors for different kinds of AGN. 
Nelson & Whittle (1996) found that there is a corre- 
lation between line width of [OIII]A5007A and stel- 
lar velocity dispersion. Wang & Lu (2001) also em- 
phasized the relation for a sample of NLSls. Re- 
cently, Greene & Ho (2005a) confirmed that the line 
width of low-ionization narrow emission lines can 
trace the stellar velocity dispersion using a sample 
of low luminous AGN selected from SDSS (York et 
al. 2000; Strauss et al. 2002; Abazajian et al. 2004). 
Zhou et al. (2006) rc-confirmcd the stronger re- 
lation for a larger sample of NLSyl selected from 
SDSS by line width of [Nil] emission line. As a test, 
Boroson (2003) have shown that the coincidence be- 
tween the BH masses estimated from the relation 
M bh — a [oiii] an d the BH masses estimated under 
the assumption of virialization for a sample of AGN 
selected from SDSS. 

As mentioned above, in order to determine the 
BH masses of the dbp emitters, we use the line width 
of [OIII] emission line as the estimator, according to 
the following equation (Gebhardt et al. 2000; Fer- 
rarese & Merritt 2001; Tremaine et al. 2002): 



M 



BH 



13±0.06^ a lOHI] 



200km • s" 



t) 



4.02±0.32 Mq (4) 



where <J[oni] is the line width of normal com- 
ponent of [OIII]A5007l. The correlation be- 
tween BH masses and EW(FeII)/EW(H/3) has 
been pointed out above as the correlation be- 
twee n( r [o/77]A5oo7i andEW(FcII)/EW(H/3) (Fig. 5) 

The anti-correlation of EW(FeII)/EW(H/3) and line 
width of [OIIIJA5007A indicates that the dbp emit- 
ters with lower BH masses have stronger Fell emis- 
sions. This result is in agreement with that of El 
analysis of single peaked BRL AGN (e. g. Sulentic, 
Marziani & Dultzin-Hacyan 2000). 

Here we compare with the correlation be- 
tween the dimcnsionlcss accretion rate and 
EW(FeII)/EW(H/3) shown in Fig 9. The di- 
mensionless ^ accretion rate rh is calculated as 

™ = m B hxi 5 38°xVs (Wandel et al. 1999; Kaspi 
et al. 2000). The Spearman Rank correlation 
coefficient is 0.23 with P nu ii ~ 24%. Such a weak 
correlation indicates that dimensionless accretion 
rate has weaker effects on Fell emission of dbp 
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emitters than BH mass. 

The mean BH masses for the 27 dbp emitters 
is 10 8 ' 05±011 Mq, the mean dimensionless accretion 
rate is io _104±011 . According to the dimensionless 
accretion rate, the accretion mode in the central re- 
gion is standard accretion rather than ADAF mode, 
because the upper limit accretion rate for ADAF 
mode is about m ~ 0.28 x a 2 where a is viscous 
coefficient a ~ 0.1 — 0.3 (Mahadevan & Quataert, 
1997; Lasota et al., 1996; Narayan et al., 1995). If 
we select a ~ 0.1 (Yi, 1996; Choi et al, 2001), all 
the dbp emitters in our sample have larger accretion 
rate than the critical value 0.0028. Furthermore, the 
27 dbp emitters have normal quasar spectra, which 
have normal big blue bumps. Thus, even if ADAF 
accretion mode dominates the accretion flow in the 
inner region of the accretion disk, especially for the 
two objects with lower accretion rate, the standard 
accretion mode should mainly dominate the bulk of 
accretion. 

4. DISCUSSIONS AND CONCLUSIONS 

We have accurately reproduced the line profiles 
in the wavelength range from 4100A to 5600A by 
the line profile of double-peaked broad Ha. The 
best fitted results indicate that the broad optical 
Fell emission lines are also double peaked and orig- 
inate from the same place where the double-peaked 
broad Balmer emission lines. More and more evi- 
dence confirms that the double-peaked broad low- 
ionization lines come from the accretion disk rather 
than an outflow or other models. In our sample, 
the double-peaked broad Ha of all the 27 dbp emit- 
ters can be best fitted by the elliptical accretion disk 
model (Eracleous et al. 1995). The best fitted re- 
sults for optical Fell emission lines indicate that Fell 
emission lines also originate from the accretion disk 
which provides the needed high density for Fell emis- 
sion lines. 

The accretion rate and the normal quasar spec- 
tra indicate that ADAF accretion mode in the dbp 
emitters with Fell emission lines is not the main ac- 
cretion mode. BH masses have more influence on the 
Fell properties of dbp emitters than dimensionless 
accretion rate. The reliable anti-correlation between 
EW(FeII)/EW(H/3) and a Ha indicate that for dbp 
emitters, the emission region of stronger Fell lines 
are far from the center black hole. We have shown 
that even for dbp emitters, the influence of orienta- 
tion on Fell intensity cannot be deconvolved from 
the other determining physical parameters of these 
AGN: mainly BH mass and, to a lesser extent, ac- 
cretion rate. 



We also get the same strong correlation between 
^5iooi anc ^ ^ Ha ( mcru des tne broad and narrow 
components near Ha) for dbp emitters in our am- 
ple as that obeyed by normal AGN, which indi- 
cates that the connection between the continuum 
and BLRs is the same for normal AGN and for dbp 
emitters. We know the luminosity of broad emis- 
sion lines, such as broad Ha, can be estimated by 
Lua = huHanennaHa^V , where, an a is the recom- 
bination coefficient for Ha, eV is the volume filled by 
gases in BLRs. The electron density in the accretion 
disk is larger than that in the normal BLRs, which 
indicates that the value of n e njyajy Q eV is a constant 
value for AGN, normal AGN and dbp emitters 

Finally we want to stress that the correlations, or 
more precisely perhaps, the trends pointed out above 
remain to be confirmed, because of two main reasons: 
First, the comparative numbers of well defined dbp 
emitters is small (which is intriguing enough), and 
second, because one of our selection criteria (the ap- 
parent strehgth of optical Fell) may bias the trends. 
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Fig. 1. The spectra of dbp emitters and the continuum under H/3 are shown in the figure. The best fitted results for 
narrow emission lines near Ha are shown in the upper right panels in each plot. The thin solid line is the spectra after 
the subtraction of the continuum. The thick solid line represents the best fitted results, the dotted line presents the 
components for [Nil] A6548, 6583A+Ha narro ™. The dash-dotted line represents the components for broad Ha. Here, 
only the first eight objects are shown in the figure. 
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Fig. 2. The best fitted results for the spectra near H/3 of dbp emitters. The dotted line represents the line spectra after 
the subtraction of the continuum. The thick solid line represents the best fitted results. The dashed line represents the 
narrow components for [OIII] doublet and/or H/3 and/or H7. The dot-dashed line represents the broad components 
for H/3, H7 and/or HeIIA4686A The double dot-dashed line represents the components for Fell emission lines in the 
wavelength range from 4200A to 5600A. Here, only the first eight objects are shown in the figure. 
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Fig. 3. The best fitted results for double-peaked broad Ha. The lower panel in each plot shows the spectrum after the 
subtraction of the broad Ha coming from the accretion disk. The model parameters are listed in Table 1. Here, only 
the first eight objects are shown in the figure. 
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Fig. 4. The correlation between EW(FeII)/EW(H/3) and 
the line width of broad Ha. The five-point star represents 
the object SDSS J2125-0813. 



Fig. 6. The correlation between EW(Fell) and the height 
ratio of [OIII]A5007A to that of H/3. The five-point star 
represents the object SDSS J2125-0813. 




Fig. 5. The correlation between EW(FeII)/EW(H/3) and 
the line width of [OIIIJA5007A. The five-point star rep- 
resents the object SDSS J2125-0813. 



Fig. 7. The correlation between the continuum luminos- 
ity and the luminosity of Ha. 
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Fig. 8. The correlation between EW(FeII)/EW(H/3) and 
the inclination angle of the accretion disk. The five-point 
star represents the object SDSS J2125-0813. 
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Fig. 9. The correlation between accretion rate and 
EW(FeII)/EW(H/3). The five-point star represents the 
object SDSS J2125-0813. 
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,98 


0. 


60u 


103 


.85+4.431 


2252 


.58 


43, 


,73 


6.98 


-0 


.43 


21 


J095744. 


,43+075124.8 





,1406 


17, 


,61 





,47u 


134 


.84+5.772 


2594 


.93 


43. 


.81 


7.44 


-0 


,81 


22 


J140502. 


,65+470747.5 





.1521 


17, 


,54 





,42u 


130 


,14+5.365 


2782 


.77 


43. 


,92 


7.37 


-0 


,63 


23 


J151535. 


,25+480530.5 





.3115 


16, 


,52 


0. 


01u 


192 


.45+6.467 


2980 


.23 


44, 


,89 


8.06 


-0 


.35 


24 


J144207. 


,69+384411.3 





.1457 


17, 


,48 


0, 


,40u 


116 


.35+2.027 


2643 


.79 


43. 


,91 


7.18 


-0 


,45 


25 


J153415. 


,41+303435.4 





.0938 


16, 


,98 


0, 


,20u 


183 


.73+2.272 


2940 


.21 


43, 


,72 


7.98 


-1 


.44 


26 


J082930. 


,59+081238.0 





,1291 


17, 


,21 


0, 


,32 


165 


.72+1.981 


3005 


.78 


43, 


,91 


7.80 


-1 


.07 



Notes: 



The first Column lists the ID number of the object. The second and third Columns present the name in the 
format of "SDSS Jhhmmss.ss+ddmmss.s" and redshift of each object. The forth column presents the PSF 
magnitude in r band of SDSS. The fifth column gives the value of R r (radio loudness) as defined by Ivezic et 
al. (2002), V represents the value of an upper limit one. The following two columns list the line width of 
[NII]A6583A and broad Ha in units of km • s _1 . The last three columns contain the value of continuum 
luminosity at 5 100 A (corrected for Galactic extinction including the redden correction, following Schlegel et 
al., 1998), the value of BH masses estimated from the line width of [OIII]A5007Jl (as a tracer of velocity 

dispersion) and the dimensionless accretion rate of the dbp emitters. 
The line width of broad Ha of J2 125-08 13 is estimated from the lincwidth of broad H/3 produced by the 
accretion disk model. Here, we just list the PSF magnitude at r band for each object. 
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TABLE 2 

PARAMETERS OF THE ACCRETION DISK MODEL 



id 


Tin 

R g 


Tout 

R g 


i 


e 






Q 






a 

km • s _1 




k 


Fell 





991 


4054 


33 





.25 




2 


.57 




2032±754 


0. 


,39 





,35 


1 


391±62 


4085±151 


89±6 





.26±0 


.01 





.98±0 


.11 


738±90 





.40 


0, 


,28 


2 


385 


969 


28 





,57 




3 


.61 




2059±400 


0, 


,43 


0, 


.13 


3 


367±37 


11946±1077 


50±3 





.24±0 


.01 


1 


,71±0 


.01 


296±25 





,37 


0. 


,11 


4 


37±17 


2660±158 


26±1 





,21±0 


.01 


1 


.29±0 


.02 


700±43 





.39 


0, 


.18 


5 


526±89 


5557±833 


54±5 





,21±0 


.01 


1 


.29±0 


.08 


771±60 


0, 


.47 


0, 


.21 


6 


526±21 


3350±255 


89 





.96±0 


.01 


2 


.24±0 


.05 


703±48 


0, 


.36 


0, 


.41 


7 


103±3 


1804±29 


21±0 





.27±0 


.02 


1 


,21±0 


.01 


915±21 


0. 


.36 


0. 


.10 


8 


2797 


i32768 


89 





.17 




7. 


,40±1 


.86 


1622±275 


0. 


,54 


0, 


.27 


9 


210±112 


5705 


47 





.03 




1 


.43±0 


.20 


1813±1078 


0. 


,34 


0, 


.42 


10 


394±40 


4953±142 


89 





.24±0 


.01 





.97±0 


.07 


1238±55 





,34 


0, 


.35 


11 


31±6 


339±40 


15±1 





.01 




1 


.58±0 


.19 


2610±270 


1. 


,00 


0. 


,53 


12 


222±18 


3382±374 


85±2 





.90±0 


.01 


2 


.03±0 


.04 


624.6±46.5 


0. 


,45 





,44 


13 


30±8 


12790±507 


25±0 





,33±0 


.01 


1 


.56±0 


.01 


292±13 


0. 


,42 


0, 


,38 


14 


408±11 


3100±116 


86±1 





.94±0 


.01 


2 


,15±0 


.02 


562.5±40.5 





,51 





,42 


15 


578±340 


^32768 


72 





.26±0 


.08 


1 


.36±0 


.01 


479±16 


0, 


,36 


0, 


,47 


16 


208±12 


4830±333 


85±2 





,84±0 


.01 


1 


.85±0 


.02 


499±24 





,44 


0, 


,36 


17 


11±1 


7779±428 


30±0 





.23±0 


.01 


1 


.44±0 


.01 


476.4±33.30 





,54 





,65 


18 


23±13 


6843±387 


23±0 





,34±0 


.01 


1 


.54±0 


.03 


461.4±27.00 


0, 


,53 


0, 


,40 


19 


88±7 


3401±226 


89±2 





.84±0 


.01 


1 


.46±0 


.02 


572.7±21.30 


0. 


.41 





.52 


20 


406±30 


4875±346 


84±7 





.87±0 


.01 


1 


.92±0 


.04 


451.5±32.40 


0. 


,51 





,53 


21 


141±43 


26691±4774 


54±7 





.28±0 


.01 


1 


.50±0 


.01 


429.6±23.70 


0, 


.30 


0, 


,23 


22 


99±19 


15609±1578 


44±2 





.32±0 


.01 


1 


.63±0 


.01 


303.9±23.70 





,51 





,75 


23 


104±16 


13542±1252 


41±2 





.33±0 


.01 


1 


.73±0 


.01 


461.7±27.30 





,54 





,42 


24 


99±13 


2551±278 


78±1 





,81±0 


.01 


1 


.54±0 


.03 


566.1±22.80 


0. 


,50 


0, 


.40 


25 


307±24 


21499±1118 


81±8 





.28±0 


.01 


1 


.56±0 


.01 


352.2±22.50 


0. 


,44 





,42 


26 


77±6 


2285±150 


76±8 





.83±0 


.01 


1 


,71±0 


.02 


676.8±29.10 


0. 


,39 


0, 


,36 


Notes: 



the error of the parameter is larger than the value of the parameter, we do not show the error. The 

number represents the name in the Table 1. 



